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SUMMARY 


Teleseismic  and  regional  data  from  nearly  100  explosions  from  the  Shagati  River  and 
Novaya  Zemlya  (USSR)  test  sites  arc  examined  to  determine  and  understand  the  spectral  nulls 
due  to  cancellation  by  pP  and  other  secondary  arrivals  and  the  scaling  of  P-wavc  spectra.  The 
U.K.  arrays  (EKA,  GBA,  and  YKa),  the  short-aperture  array  NORESS  and  stations  MAJO 
and  NORSAR  (center  element)  provided  the  teleseismic  data,  whereas  the  regional  data  came 
from  the  CDSN  station,  WMQ.  Results  from  the  array  data,  derived  by  a  least  squares  inver¬ 
sion  that  isolates  the  source  and  receiver  terms  for  each  frequency,  should  be  relatively  free 
from  frequency-dependent  station  effects.  A  combination  of  several  techniques,  including 
spectral  ratios  P/P-coda,  Pn/Lg,  inter-shot  ratios  of  P  and  Pn,  and  multichannel  deconvolution, 
provide  effective  means  for  determination  and  interpretation  of  spectral  nulls.  Analysis  of 
data  from  nearly  80  Shagan  River  explosions  suggests  that,  in  general,  the  spectra  of  P  waves 
have  three  different  types  of  spectral  nulls:  (1)  well-defined  nulls  due  to  cancel  1  -non  by  pP  at 
the  expected  frequency  of  about  4  Hz  in  both  teleseismic  aid  regional  data  to;  1  trge  (mb  =  6) 
shots,  (2)  low-frequency  (less  than  about  2  Hz)  nulls  m  tc’eseismic  data  tfru  appear  to  be 
related  to  both  local  geology  and  depth  of  burial,  and  (3)  distinct  nulls  at  about  3  Hz  in 
teleseismic  data  for  large  and  small  shots,  probably  due  to  ••deseismic  pat  >  effect  such  as  a 
multiple  arrival.  A  comparison  of  P  and  P-coda  shows  that  the  low-fi  equency  P-coda  is 
mainly  due  to  the  scattering  of  explosion-genera'ed  hg  and  the  spectral  rado  P/P-coda  may  be 
used  to  estimate  shot  depths.  Scaling  of  P-wave  spectra  is  investigated  over  the  frequency 
range  of  0.5  to  10  Hz  by  comparison  of  observed  data  from  WMQ  and  NORSAR  with  syn¬ 
thetics  assuming  several  t*  and  pP  reflection  coefficient  values.  In  qualitative  agreement  with 
Stewart’s  (1988)  analysis  of  Shagan  River  explosions,  the  observed  P  spectra  vary  more 
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slowly  with  mb  than  predicted  by  either  Mueller  and  Murphy  (1971)  or  von  Seggem  and 
Blandford  (1972)  scaling  relationships  and  the  differences  are  more  pronounced  at  lower  fre¬ 
quencies  than  at  the  higher  frequencies. 


Teleseismic  data  from  19  Novaya  Zemlya  explosions  recorded  at  the  EKA,  GBA,  and 
YKA  arrays  are  analyzed  to  understand  the  characteristics  of  the  pP  arrival.  Large  (mb  ~  6) 
shots  indicate  distinct  spectral  nulls  at  about  2.5-3.0  Hz,  as  expected  for  nulls  due  to  pP. 
Spectral  ratio  of  P/P-coda,  averaged  over  all  available  sensors  at  all  arrays,  and  inter-shot 
spectral  ratios  of  large  and  small  shots  at  common  receivers,  provide  reliable  methods  for 
estimating  spectral  nulls.  A  comparison  of  synthetics  with  the  observed  data  suggests  the 
amplitude  ratio  pP/P  to  be  about  0.5  for  large  shots.  These  results  obtained  by  using  new 
frequency-domain  techniques  are  helpful  in  resolving  the  controversy  over  the  pP  arrival  time 
and  amplitude. 


IV 


TABLE  OF  CONTENTS 


Page 


1.  SUMMARY  iii 

2.  INTRODUCTION  1 

3.  STUDY  OF  SHAGAN  RIVER  EXPLOSIONS  5 

3.1.  ANALYSIS  OF  EK A  ARRAY  DATA  5 

3.1.1.  Spectral  Nulls  in  Source  Spectra  S 

3.1.2.  Low-Frequency  Nulls  and  Near-Source  Geology  9 

3.1.3.  Comparison  of  P  and  P  Coda  and  Origin  of  Secondary  Arrivals  16 

3.2.  ANALYSIS  OF  NORESS  ARRAY  DATA  22 

3.3.  ANALYSIS  OF  MULTI-STATION  DATA  FOR  USSR  JVE  SHOT  23 

3.4.  ANALYSIS  OF  REGIONAL  DATA  AT  WMQ  26 

3.5.  SCALING  OF  P- WAVE  SPECTRA  34 


3.5.1.  Comparison  of  Theory  and  Observed  Regional  Data 

3.5.2.  Comparison  of  Theory  and  Observed  Teleseismic  Data 

3.6.  DISCUSSION 

4.  STUDY  OF  NOVAYA  ZEMLYA  EXPLOSIONS 

4.1.  P-WAVE  SOURCE  SPECTRA  FROM  U.K.  ARRAY  DATA 

4.2.  INTER-SHOT  RATIOS  AND  ESTIMATES  OF  NULL  FREQUENCY 

4.3.  ESTIMATES  OF  pP/P  BY  COMPARISON  WITH  SYNTHETICS 

5.  CONCLUSIONS 

6.  ACKNOWLEDGMENTS 

7.  REFERENCES 


Accesion  For 


NTIS  CRA&I 
DTIC  TAB  □ 

U..a.vounced  i  ; 
Justification 

By  . 

Out.  ibutioii/ 

Availability  Cv-.o 


Aval.  u.  ..  i  or 
D'st  !  S:\  -ial 


M 


34 

37 

45 

51 

51 

58 

66 

71 

73 

74 


2.  INTRODUCTION 


An  understanding  of  the  spectral  nulls  in  the  P-wave  spectra  of  underground  nuclear 
explosions  is  essential  for  constraining  the  source  depth  and  for  determining  the  bias  in  the 
body  wave  magnitude  resulting  from  the  constructive  or  destructive  interference  with  pP  and 
later  arrivals.  Determination  of  the  correct  spectral  nulls  due  to  pP  for  Shagan  River  explo¬ 
sions  is  especially  important  because  of  the  complete  lack  of  agreement  among  various  work¬ 
ers.  The  multichannel  deconvolution  of  array  data  by  Der  etal.  (1987a)  provided  pP-P  or 
delay  time  values  of  0.55  sec  for  Kazakh  events  recorded  at  NORSAR  and  0.4  sec  for  those 
at  EKA.  Using  network-averaged  P-wave  spectra,  Murphy  et  al.  (1991)  obtained  delay  time 
of  0.82  sec  for  an  explosion  with  yield  of  119  kt  at  Shagan  River  test  site.  Stewart  (1988) 
examined  62  Shagan  River  explosions  recorded  at  the  four  U.K.  arrays  (EKA,  GBA,  WRA, 
and  YKA),  and  the  deconvolved  wavefomts  suggest  that  pP  delay  times  for  large  (mb  =  6) 
explosions  vary  between  about  0.4  to  0.9  sec.  All  these  values  of  delay  times  are  inconsistent 
with  what  is  known  about  the  site  conditions  and  testing  practices  which  would  suggest  a 
delay  time  for  large  explosions  to  be  no  larger  than  about  0.25  sec.  An  examination  of  the 
pulse  durations  and  rise  times  of  Shagan  River  explosions  by  Stewart  (1988)  indicated  that  the 
source  corner  frequency  scale-.,  at  most,  as  Y~l/l°,  which  is  much  smaller  than  the  Y_,/3 
predicted  by  most  source  models  and  observed  by  Lyman  etal.  (1986)  in  the  rise  times  of 
NTS  explosions.  The  main  objectives  of  this  study  are  to  (1)  determine  and  understand  spec¬ 
tral  nulls  due  to  cancellation  by  pP  am  other  secondary  arrivals  in  P-wave  spectra  of  Shagan 
River  (USSR)  explosions,  and  (2)  compare  the  observed  variation  of  P-wave  spectra  as  a 
function  of  mb  with  von  Seggern  and  Blandford  (1972)  and  Mueller  and  Murphy  (1971)  scal¬ 
ing  relationships.  We  examined  teleseismic  and  regional  data  from  79  underground  nuclear 
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explosions  at  the  Shagan  River  test  site  (Table*  1).  The  list  includes  40  events  (28  from  SW 
and  12  from  NE  region)  recorded  at  the  EKA  array,  33  shots  at  MAJO,  25  at  NORSAR  (sin¬ 
gle  channel),  20  at  the  NORESS  array,  and  17  at  WMQ.  Most  events  are  recorded  at  more 
than  one  a’-ray  or  station. 

A  study  of  characteristics  of  the  P-wave  spectra  of  Novaya  Zemlya  explosions  is  impor¬ 
tant  because  of  the  possibility  of  all  future  USSR  shots  taking  place  in  this  region.  More¬ 
over,  there  is  significant  discrepancy  between  the  results  of  alternate  approaches  to  under¬ 
standing  the  character  of  the  pP  reflection  from  these  shots.  As  pointed  out  by  Burdick 
(1990),  the  arrival  time  of  "effective  pP"  is  substantially  later  than  one  would  predict  for  a 
point  source  within  an  elastic  medium  and  estimates  of  its  amplitude  differ  by  as  much  as  a 
factor  of  5.  Our  analysis  of  digital  data  for  19  Novaya  Zemlya  explosions  recorded  at  the 
U.K.  arrays  provides  information  regarding  pP  that  appears  to  be  consistent  with  local  geology 
of  the  test  site. 
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TABLE  1 


79  SHAG  AN  RIVER  EXPLOSIONS  USED  IN  STUDY 


No. 

DATE 

LAT* 

LON* 

mb* 

ERA 

MAJO 

NA0 

1 

15  Jan  1965 

49.940 

79.010 

5.87 

X 

_ 

2 

19  Jun  1968 

49.982 

79.003 

5.28 

X 

- 

3 

30  Nov  1969 

49.913 

78.961 

6.02 

X 

• 

4 

30  Jun  1971 

49.949 

78.986 

4.94 

X 

X 

5 

31  May  1974 

49.950 

78.852 

5.81 

X 

- 

6 

16  Oct  1974 

49.979 

78.898 

5.41 

X 

- 

- 

7 

27  Dec  1974 

49.943 

79.011 

5.50 

X 

X 

8 

27  Apr  1975 

49.949 

78.926 

5.51 

X 

- 

9 

30  Jun  1975 

50.004 

78.957 

4.52 

X 

X 

10 

29  Oct  1975 

49.946 

78.878 

5.61 

X 

. 

11 

25  Dec  1975 

50.044 

78.814 

5.69 

X 

. 

12 

21  Apr  1976 

49.890 

78.827 

5.12 

• 

X 

13 

09  Jun  1976 

49.989 

79.022 

5.07 

- 

X 

14 

04  Jul  1976 

49.909 

78.91 1 

5.85 

X 

- 

15 

28  Aug  1976 

49.969 

78.930 

5.74 

X 

- 

16 

07  Dee  1976 

49.922 

78.816 

5.80 

X 

X 

!7 

29  May  1977 

49.937 

78.7  0 

5.75 

X 

X 

18 

29  Jun  1977 

50.006 

78.8  .9 

5.20 

X 

X 

19 

05  Sep  1977 

50.035 

78.9  1 

5.73 

X 

X 

20 

30  Nov  1977 

49.958 

78.8  i5 

5.89 

X 

• 

21 

11  Jun  1978 

49.898 

78.7  >7 

5.83 

X 

X 

X 

22 

05  Jul  1978 

49.887 

78.8  1 

5.77 

X 

X 

X 

23 

15  Sep  1978 

49.916 

78.8  9 

5.89 

X 

X 

X 

24 

04  Nov  1978 

50.034 

78.943 

5.56 

X 

X 

X 

25 

23  Jun  1979 

49.903 

78.855 

6.16 

X 

X 

X 

26 

07  Jul  1979 

50.026 

78.991 

5.84 

X 

X 

X 

27 

04  Aug  1979 

49.8(<4 

78.904 

6.13 

X 

X 

X 

28 

18  Aug  1979 

49.9  »3 

78.938 

6.13 

X 

X 

• 

29 

28  Oct  1979 

49.9  '3 

78.997 

5.98 

X 

X 

X 

30 

02  Dec  1979 

49.891 

78.796 

5.99 

X 

X 

X 

31 

23  Dee  1979 

49.916 

78.755 

6.13 

X 

X 

- 

32 

12  Jun  1980 

49.98 ) 

79.001 

5.52 

- 

X 

- 

33 

29  Jun  1980 

49.93  1 

78.815 

5.69 

X 

X 

- 

34 

14  Sep  1980 

49.92 . 

78.802 

6.21 

- 

X 

X 

35 

12  Oct  1980 

49.961 

79.028 

5.88 

X 

X 

X 

36 

14  Dee  1980 

49.899 

78.9.18 

5.93 

X 

X 

X 

37 

27  Dec  1980 

50.057 

78.98 1 

5.87 

X 

X 

X 

38 

29  Mar  1981 

50.007 

78.98? 

5.49 

X 

X 

X 

39 

22  Apr  1981 

49.885 

78.81 ) 

5.94 

X 

X 

X 

40 

27  May  1981 

49.985 

78.98 ) 

5.30 

- 

- 

X 

41 

13  Sep  1981 

49.910 

78.91  > 

6.06 

X 

X 

- 

NRSS 


WMQ 
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42 

18  Oct  1981 

49.923 

78.859 

6.00 

X 

X 

43 

29  Nov  1981 

49.887 

78.860 

5.62 

X 

m 

44 

27  Dec  1981 

49.923 

78.795 

6.16 

• 

X 

45 

25  Apr  1982 

49.903 

78.913 

6.03 

• 

X 

46 

31  Aug  1982 

49.924 

78.761 

5.20 

X 

47 

26  Dee  1982 

50.071 

78.988 

5.58 

X 

„ 

48 

12  Jun  1983 

49.913 

78.916 

6.02 

, 

X 

49 

26  Oct  1983 

49.901 

78.828 

6.04 

• 

X 

50 

25  Apr  1984 

49.929 

78.870 

5.90 

. 

X 

51 

26  May  1984 

49.969 

79.006 

6.01 

X 

52 

14  Jul  1984 

49.893 

78.884 

6.10 

X 

53 

27  Oct  1984 

49.920 

78.777 

6.19 

X 

54 

02  Dec  1984 

49.989 

79.011 

5.77 

X 

55 

16  Dec  1984 

49.926 

78.820 

6.12 

X 

56 

28  Dee  1984 

49.866 

78.703 

6.00 

X 

57 

10  Feb  1985 

49.888 

78.781 

5.83 

X 

58 

25  Apr  1985 

49.914 

78.902 

5.84 

59 

15  Jun  1985 

49.898 

78.845 

6.05 

60 

30  Jun  1985 

49.848 

78.658 

5.92 

61 

20  Jul  1985 

49.936 

78.785 

5.89 

62 

12  Mar  1987 

49.939 

78.823 

5.3) 

63 

03  Apr  1987 

49.928 

78.829 

6.12 

64 

17  Apr  1987 

49.886 

78.691 

5.92 

65 

20  Jun  1987 

49.913 

78.735 

6.03 

66 

02  Aug  1987 

49.880 

78.917 

5.83 

67 

15  Nov  1987 

49.871 

78.791 

5.98 

68 

13  Dec  1987 

49.989 

78.844 

6.06 

69 

27  Da*  1987 

49.864 

78.758 

6.00 

70 

13  Rb  1988 

49.954 

78.910 

5.97 

71 

03  Apr  1988 

49.917 

78.945 

5.99 

72 

04  May  1988 

49.928 

78.769 

6.09 

73 

14  Jun  1988 

50.045 

79.005 

4.80 

74 

14  Sep  1988 

49.870 

78.820 

6.03 

75 

12  Nov  1988 

50.056 

78.991 

5.20 

76 

17  Dee  1988 

49.818 

78.910 

5.80 

77 

12  Feb  1989 

49.930 

78.740 

5.90 

78 

08  Jul  1989 

49.870 

78.820 

5.60 

. 

79 

02  Sep  1989 

50.020 

79.050 

500 

X 

X 

X 


*  from  Ringdal  and  Marshall  (1989),  Norsar  Repi.  2-88/89 
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3.  STUDY  OF  SHAGAN  RIVER  EXPLOSIONS 


3.1.  ANALYSIS  OF  EKA  ARRAY  DATA 

We  examined  the  short-period,  vertical-component  data  from  40  Shagan  River  explosions 
well  recorded  at  the  Eskdalemuir  (EKA)  array,  which  has  20  elements  arranged  in  an  "L" 
configuration  with  element  spacing  of  0.9  km  and  maximum  spacing  of  9.8  km  (Bache  et  at. , 
1985).  The  epicentral  distance  of  Shagan  explosions  to  the  array  is  about  47°  along  the 
azimuth  of  about  310°  or  N50°W.  As  compared  to  results  from  a  single  receiver,  use  of  mul¬ 
tichannel  data  has  the  advantage  that  frequency-dependent  receiver  effects,  which  can  drasti¬ 
cally  influence  the  observed  signals  (e.  g. ,  Murphy  et  at. ,  1989),  can  be  eliminated  or  at  least 
suppressed. 

3.1.1.  Spectral  Nulls  In  Source  Spectra 

Using  a  window  length  of  6.4  sec.  starting  2  sec  before  the  onset  of  direct  P  and  with 
Parzen  taper,  Fourier  spectra  were  obtained  for  all  data  except  a  few  channels  with  spikes  or 
clipping.  These  P  windows  are  effectively  only  a  few  sec  long  and  are  kept  short  so  as  to 
minimize  the  effects  of  later  secondary  arrivals.  The  spectra  were  corrected  for  instrument 
response  and  for  noise  by  taking  Fourier  spectra  of  a  sample  of  noise  preceding  the  P  arrival. 
A  least  squares  inversion  method  that  isolates  the  source  and  receiver  terms  for  each  fre¬ 
quency  was  applied  to  all  spectral  amplitudes  with  the  signal/noise  (or  S/N)  power  ratio 
greater  than  1.  The  derived  source  terms,  arranged  >n  order  of  decreasing  m^  for  28  shots 
from  the  southwest  Shagan  test  site  region,  are  shown  in  Figures  la,  b  whereas  the  remaining 
12  from  the  northeast  Shagan  region  are  shown  in  Figure  1c.  All  (40)  source  spectra, 
arranged  in  the  same  order  as  in  Figure  1  but  with  incremental  shifts  along  both  axes,  are 
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FREQUENCY  (HZ)  FREQUENCY  (HZ) 


Figure  la.  Source  spectra,  based  on  the  P(6.4  sec)  window,  collected  for  instrument  response, 
of  the  14  largest  SW  Shagan  shots  recorded  at  EKA,  arranged  in  order  of  decreasing  mb. 
Mean  values  of  the  amplitude  ratio  P/P-coda  (0.5-2.0  Hz)  are  also  indicated.  Prominent  low- 
frequency  nulls  are  generally  associated  with  smaller  values  of  P/P-coda. 
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LOG  AMPLITUDE 

1.00  2.00  3.00  4.00  5.00  6.00  7.00  8.00 


FREQUENCY  (HZ)  FREQUENCY  (HZ) 

Figure  lb.  Similar  to  Figure  la  for  the  next  14  SW  Shagan  shots,  arranged  in  order  of 
decreasing  mb.  The  low-frequency  n.  Us,  more  distinct  for  smaller  explosions,  are  generally 
associated  with  smaller  values  of  P/P-c  xia. 
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Figure  lc.  Similar  to  Figure  la  for  12  NE  Shagan  shots,  arranged  in  order  of  decreasing  mb. 
The  low-frequency  nulls,  generally  more  distinct  for  smaller  explosions  and  for  shots  in  the 
NE  Shagan  region,  are  generally  associated  with  smaller  values  of  P/P-coda. 
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shown  in  Figure  2.  An  examination  of  the  spectra  in  Figures  1  and  2  suggests  three  general 
trends:  (1)  a  majority  of  shots,  including  the  cratering  shot  of  15  January  1965,  have  spectral 
nulls  at  about  3  Hz,  (2)  several  explosions  have  spectral  nulls  at  about  4  Hz,  and  (3)  a  major¬ 
ity  of  shots  have  low-frequency  (less  than  about  2  Hz)  nulls  which  are  generally  more  distinct 
for  smaller  explosions  and  for  shots  in  the  northeast  Shagan  region.  These  three  types  of 
spectral  nulls  are  also  evident  in  the  average  po’^er  spectrum  for  all  40  shots  (Figure  3)  which 
indicates  the  most  pronounced  null  to  be  at  about  3  Hz  and  less  distinct  nulls  at  about  4  Hz 
and  1.5  Hz. 

3.1.2.  Low-Frequency  Nulls  and  Near-Source  Geology 

The  Shagan  River  test  site  is  believed  to  be  underlain  by  abou  100  m  of  unconsolidated 
sediments  with  the  thic  kness  decreasing  gradually  towards  southwest  (Bonham  et  al. ,  1980). 
The  P(6.4  sec)  source  spectra  of  10  explosions  (a  subset  of  the  40  source  spectra  in  Figure  I 
but  without  correction  for  instrument  response)  with  shot  points  lying  approximately  along  a 
SW-NE  profile  across  the  test  site  (Figure  4a)  suggests  a  gradual  build-up  of  the  low- 
frequency  nulls  as  the  thickness  of  sediments  increases.  Moreover,  the  peak-amplitude  and 
the  null  frequencies  appear  to  be  decreasing  with  increasing  thickness  of  sediments.  A  null  at 
about  3  Hz  is  also  observed  in  most  spectra.  Figures  la  and  lb  indicate  that  these  low- 
frequency  nulls  are  als  >  associated  with  smaller  mb  and  therefore  shallower  shots  in  the  SW 
Shagan  region.  Source  spectra  derived  from  the  P(3.2  sec)  windows  for  the  same  10  explo¬ 
sions,  shown  in  Figure  4b,  indicate  almost  total  absence  of  the  low-frequency  nulls  and  some 
evidence  of  nulls  at  about  4  and  3  Hz.  Moreover,  ft  r  most  1  Jge  shots,  the  4  Hz  null  is  rela¬ 
tively  more  prominent  in  spectra  of  the  shorter  window  (Figure  4b)  than  in  spectra  of  the 
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Figure  2.  Source  spectra  of  40  explosions,  same  as  in  Figures  la,  b,  c  and  arranged  in  the 
same  order  but  with  incremental  shifts  along  both  axes.  The  log  amplitude  values  are  normal¬ 
ized  to  the  first  frequency  and  the  arrow  divides  the  SW  and  NE  Shagan  populations. 
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1  2  3  4  5  6 

FREQUENCY  (HZ) 


Figure  3.  Average  power  s  >ectrum  (log  units)  of  all  40  shots  recorded  at  EKA  indicating  a 
distinct  null  at  about  3  Hz  a  id  weaker  nulls  at  about  4  Hz  and  1.5  Hz. 
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(a)  P(6.4  sec) 


(b)  P(3.2  sec) 


FREQUENCY  (HZ)  FREQUENC  Y  (HZ) 

Figure  4.  Source  spectra  derived  from  EKA  array  data,  not  corrected  for  instrument  response, 
of  10  explosions  with  shot  points  lying  approximately  along  a  SW-NE  profile  across  the  test 
sire,  for  (a)  P(6.4  sec)  and  (b)  P(3.2  sec)  windows.  Note  the  gradual  build-up  of  the  low- 
frequency  nulls  which  are  more  dominant  in  (a)  than  in  (b).  Numbers  (in  parentheses)  next  to 
the  spectra  indicate  mb. 
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longer  window  (Figure  4a).  Even  the  3  Hz  null  is  relatively  less  prominent  in  Figure  4b  than 
in  Figure  4a,  suggesting  the  possibility  that  its  origin  is  also  due  to  a  somewhat  later  arrival. 
It  is  therefore  clear  that  the  low-frequency  nulls  an-  due  to  secondary  (later)  arrivals.  Further¬ 
more,  a  probable  reason  for  their  existence  is  gr  ater  near-source  scattering  for  explosions 
under  thicker  low  velocity  sediments  or  shallower  shot  depths  because  the  near-source 
environment  in  both  cases  is  likely  to  be  associated  with  greater  near-source  heterogeneity. 
Thin  low  velocity  sediments  may  not  be  "seen"  by  seismic  waves  of  low  frequency  or  long 
wavelengths  and  therefore  not  make  significant  contribution  to  the  scattered  waves. 

Figure  5  shows  source  deconvolutions  for  a  set  of  Shagan  River  explosions  recorded  at 
EKA  as  d  rived  by  Der  et  al.  (1987b).  The  first  three  shots  are  from  the  SW  and  the  remain¬ 
ing  three  from  the  NE  region  of  the  test  site.  The  corresponding  source  spectra,  a  subset  of 
those  in  Figure  1  but  without  the  correction  for  instrument  response  (similar  to  Figure  4a),  are 
shown  in  Figure  6.  Similar  to  the  results  in  Figure  4,  a  comparison  of  Figures  5  and  6  also 
suggests  that  the  low-frequency  nulls  are  associated  with  prominent  secondary  arrivals  within 
a  few  sec  of  the  first  P.  The  deconvolutions  have  distinct  low-frequency  secondary  arrivals 
within  a  few  sec  of  the  first  P  for  the  last  4  shots  in  Figure  5  and  spectra  of  the  same  4  shots 
show  distinct  low-frequency  nulls  (Figure  6).  The  three  SW  shots  in  Figure  5  have  fairly  dis¬ 
tinct  pP  phases  and  the  spectra  of  the  first  two  show  the  expected  nulls  due  to  cancellation  by 
pP  at  about  4  Hz  (Figure  6).  The  third  (29  June  80)  shot  indicates  a  strong  secondary  arrival 
with  its  positive  peak  arriving  about  1.3  sec  after  the  positive  peak  in  the  first  P.  These  two 
same-polarity  arrivals  should  give  rise  to  spectral  peaks  at  about  0.8  Hz,  1.6  Hz,  2.4  Hz,  .  .  , 
in  good  agreement  with  the  observed  source  spectra  in  Figure  6.  Similar  interpretations 
appear  to  be  valid  for  the  strong  low-frequency  peaks  and  nulls  in  the  spectra  of  the  three  NE 
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2  DEC  79  (SW),  6.0,  0.605 


4.0  sec 


Figure  5.  Deconvolved  source  functions  of  6  Shagan  River  explosions  derived  from  EKA 
array  data.  The  P  and  pP  arrivals,  if  present,  are  denoted  by  arrows.  mb  and  mean  values  of 
the  amplitude  ratio  P/P-coda  (0.5-2.0  Hz)  are  also  indicated. 
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Figure  6.  Source  spectra,  not  corrected  for  instrument  response,  of  the  same  6  Shagan  River 
explosions  as  in  Figure  5. 
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shots  which  do.  not  show  distinct;  pP  arrivals.  Bach  set  of  3  shots  is  closely  spaced  (within 
about  5  km),  and  the  codas  from  each  set  seem  to.  have  several  common  arrivals. 


3.0.  Comparison  of  P  and  P  Coda  and  Origin  of  Secondary  Arrivals 


In  order  to  understand  the  origin  of  secondary  arrivals  and  the  associated  spectral  nulls, 
we  made  a  comparison  of  the  spectral  characteristics  of  P  and  P  coda.  It  is  known  that  the 
null  due  to  cancellation  by  pP  in  the  P-wave  spectra  is  significantly  reduced  in  the  spectra  of 
P  coda  (Gupta  and  Blandfotd,  1987).  Furthermore,  the  amplitude  ratio  P/P-coda  may  provide 
an  indication  of  near-source  scattering,  such  as  scattering  of  the  explosion-generated  Rg  into 
teleseismic  P  arrivals  (Gupta  et  al.  ,1991). 

For  an  explosion  source  with  source  function  S(f)  within  a  homogeneous  half  space,  the 
fundamental-mode  Rayleigh  wave  amplitude  Rg(f)  is  given  by  the  e  ^pression  (Hudson  and 
Douglas,  1975): 


-4.4  -  f 


Re  (0  =  f15 10  0 


r0-5  p  a3-5 


(1) 


where  <t»((o)  is  the  Fourier  transform  of  the  source  RDP  <J>(t),  h  is  the  shot  depth,  to  =  2rcf,  a  is 


the  compressional-wave  velocity  of  the  shot  medium,  p  is  the  density,  r  is  the  epicentral  dis¬ 
tance,  A  is  a  numerical  constant,  and  a  Poisson’s  ratio  of  0.25  is  assumed.  The  far- field  P- 


wave  displacement,  in  the  frequency  domain,  may  be  written  as  (Hudson  and  Douglas,  1975): 


P(f)  = 


1 


<0  <{>(co) 


(2) 


4  n  p  a3  R 

where  R  ;s  the  te'eseismic  source-receiver  distance.  If  the  source-receiver  path  and  the 


medium  velocity  do  not  change,  equations  (1)  and  (2)  yield 
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.h 

a 


(3) 


(0  =  k  C(f)  f 0,5  10~4A 

where  k  is  a  constant  and  C(f)  denotes  the  scattering  function  Rg->P.  Assuming  that  the 
low-frequency  P-coda  in  teleseismic  data  is  mainly  due  to  the  scattering  of  Rg-»P,  we  replace 
Rg-*P  and  P  in  equation  (3)  by  P-coda  and  P,  respectively,  and  obtain 

log  *^T  (0  =  4.4-  f- log  k  -  log  C(0  -  0.5  log  f  (4) 

This  shows  that  if  we  consider  explosions  for  which  the  Rg-»P-coda  scattering  function  may 
be  assumed  to  be  the  same  (such  as  closely  spaced  explosions)  and  the  frequency  is  kept  fixed 
(i.e.  for  a  narrow  frequency  band),  a  plot  of  log  P/P-coda  versus  shot  depth,  h  should  have  a 
slope  of  about  4.4  f/a.  Observations  of  he  spectra  of  Rg  signals  from  both  explosions  and 
very  shallow-focus  earthquakes  indicate  the  dominant  energy  in  Rg  to  be  confined  to  frequen¬ 
cies  less  than  about  2  Hz  and  the  most  important  parameter  for  the  generation  of  Rg  to  be 
source  depth  (e.g.  Kafka,  1990). 

Excluding  a  few  channels  with  spikes  or  clipped  data,  spectra  and  spectral  ratios  P/P- 
coda  were  obtained  for  each  rec  >rd  by  selecting  a  window  length  of  12.8  sec  (beginning  4  sec 
before  the  onset  of  P)  as  the  P-  window  and  the  following  25.6  sec  signal  as  the  P-coda,  and 
applying  Parzen  taper  to  each  (see  Figure  1  of  Gupta  and  Blandford,  1987  for  examples).  A 
correction  for  noise  was  made  by  selecting  a  sample  of  noise  before  the  onset  of  P  and 
obtaining  its  spectra.  The  average  amplitude  ratio  of  P/P-coda.  /er  the  frequency  range  of 
0.5-2.0  Hz,  was  obtained  by  using  only  those  data  points  for  which  the  S/N  power  ratio  was 
at  least  2  for  each  of  P  and  P-coda.  All  available  amplitude  ratios  P/P-coda  were  used  as 
input  to  a  least  squares  inversion  that  separated  the  source  and  receiver  terms.  A  plot  of  the 
source-term  amplitude  ratio  P/P-coda  versus  mb  (Figure  7a)  suggests  generally  larger  P/P-coda 
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Figure  7.  Source  term  amplitude  ratio  P/P-coda  (log  units),  averaged  over  the  frequency 
range  of  0.5-2.0  Hz,  for  40  Shagan  River  explosions  recorded  at  th  •  EKA  array  plotted  versus 
(a)  mb  and  (b)  estimated  shot  depth.  The  least  squ.ires  linear  regi  jssion  (dashed  line),  mean 
slope  (with  associated  standard  deviation),  and  one  standard  deviation  of  residuals  (SD)  are 
indicated  on  each  plot.  Six  pairs  of  closely  spaced  explosions  are  identified  by  the  letters  A 
through  F  and  the  SW  and  NE  explosion  populations  are  denoted  by  X  and  O,  respectively. 
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for  larger  explosions.  Figure  7a  also  shows  that  the  NE  Shagan  shots  generally  have  smaller 
values  of  P/P-coda  than  the  SW  shots.  This  will  be  consistent  with  generally  larger  P-coda 
and  greater  complexity  observed  for  the  NE  shots  than  for  the  SW  shots  by  Bache  et  al. 
(1984)  and  Marshall  et  al.  (1985). 

A  relationship  between  the  average  ratio  P/P-coda  and  shot  depth  was  obtained  by  con¬ 
verting  the  known  mb  into  shot  depth  by  using  empirical  relationships  based  on  the  informa¬ 
tion  on  shot  depths  and  yields  recently  released  by  Bocharov  et  al.  (1989).  These  relation¬ 
ships,  derived  from  the  Shagan  River  and  Konystan  (Murzhik)  explosions  are  (R.  S.  Jih.  writ¬ 
ten  communications,  1990): 

mb  =  0.690  log  Y  +  4.60.5  (5) 

where  Y  is  the  yield  in  kt,  and 

log  h  =  0.241  log  Y  + 2.174  (6) 

where  h  is  shot  depth  in  meters.  These  two  equations  would  indicate  an  explosion  of  mb  =  6 

to  have  a  yield  of  about  105  kt  and  a  depth  of  458  m.  Using  equations  (5)  and  (6),  one 

obtains 

log  h  =  0.349  mb  t-  0.567  (7) 

which  can  be  used  to  estimate  shot  depth  for  an  explosion  with  known  mb.  A  plot  of  the 

observed  amplitude  ratio  P/P-coda  versus  estimated  shot  depth  in  km  (Figure  7b)  indicates  a 

mean  slope  of  0.943  with  one  standard  deviation  of  residuals  of  only  0.095. 

Assuming  the  P-wave  velocity  for  the  uppermost  0.5  km  of  the  crust  in  the  Shagan  test 
site  region  to  be  5  km/sec,  4.4  f/a  =  1.1,  if  f  is  taken  to  be  the  average  of  0.5  and  2.0  or  1.25 
Hz.  Good  agreement  between  the  mean  slope  expected  from  simple  theory  and  observation 
(Figure  7b)  suggests  that  the  low-frequency  P-coda  is  mainly  due  to  the  scattering  of 
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explosion-generated  Rg. 

Figure  7  includes  results  from  six  (including  two  from  NE  Shagan)  pairs  of  explosions 
\vithin  less  than  2  km  of  each  other  (identified  by  the  letters  A  through  F).  It  is  known  that 
propagation  paths  can  exert  significant  influence  on  seismic  recordings  (e.g.  Blandford,  1981) 
so  that  results  from  closely  spaced  explosions  recorded  at  a  common  station  should  be  con¬ 
sidered  more  reliable  than  others.  It  is  important  to  note  that  the  slope  from  each  of  the  six 
pairs  of  shots  in  Figure  7b  is  nearly  the  same  as  derived  from  the  average  of  all  (40)  shots 
and  is  in  good  agreement  with  the  theoretically  expected  value. 

Results  similar  to  those  in  Figure  7  but  for  the  fr«  quency  range  of  3.0-5.0  Hz  are  shown 
in  Figure  8  which  shows  good  separation  between  explosions  from  the  SW  and  NE  regions  of 
the  Shagan  River  test  site.  The  lack  of  any  systematic  variation  with  shot  depth  is  consistent 
with  the  suggestion  of  explosion-generated  Rg  contributing  to  P-coda  because  Rg  is  important 
on  y  for  frequencies  less  than  about  2  Hz.  Following  Levander  and  Hill  (1985),  the  relatively 
larger  high-frequency  P  coda  for  shots  in  the  NE  Shagan  region  could  be  due  to  the  presence 
of  an  irregular  low-velocity  surface  layer  in  that  region. 

If  the  scattering  of  explosion-generated  Rg  contributes  mainly  to  P  coda,  what  is  the  ori¬ 
gin  of  the  secondary  arrivals  appearing  within  a  few  sec  of  the  first  P  and  causing  distinct 
low-frequency  nulls  in  the  spectra  of  the  initial  P  window?  The  answer  appears  to  be  that 
these  secondary  arrivals  are  due  to  near-source  or  local  scattering  because  our  observations 
suggest  the  low-frequency  nulls  to  be  stronger  for  shots  associa  ed  with  greater  near-source 
heterogeneity.  The  scattering  may  include  body-wave  (P-to-P  and  S-to-P)  scattering  and  the 
near-source  S  may  be  generated  by  mechanisms  such  as  source  anisotropy  or  strain  release. 
Local  scattering,  expected  to  be  larger  for  both  shallower  shots  and  for  thcce  in  NE  Shagan, 
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40  SR  SHOTS,  EKA  ARRAY  (3.0-5.0  Hz) 
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Figure  8.  Similar  to  Figure  7  but  for  the  frequency  range  of  3.0-5.0  Hz.  There  is  no 
significant  dependence  of  P/P-coda  on  mb  or  estimated  shot  depth  but  the  SW  and  NE  explo¬ 
sion  populations  appear  well  separated. 
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may  therefore  be  responsible  for  the  low-frequency  nulls. 

P  coda  is  known  to  be  an  excellent  measure  of )  ields  of  underground  nuclear  explosions 
at  both  NTS  and  USSR  test  sites  (Gupta  etal ,  1985).  The  observed  stability  of  P  coda 
appears  to  imply  that  the  Rg-to-P  scattering  must  be  occurring  over  a  fairly  large  region  of  the 
shallow  crust.  Our  observations  and  analysis  suggest  that  the  most  important  parameter  con¬ 
trolling  Rg-to-P  scattering  within  a  given  region  is  the  shot  depth.  Thus  local  body-wave 
scattering  may  be  responsible  for  the  low-frequency  nulls  whereas  the  scattering  of  Rg-to-P 
occurs  on  a  regional  basis  and  appears  in  P  coda.  A  prominent  low-frequency  spectral  null 
and  relatively  large  value  of  P/P-coda  observed  for  the  nuclear  shot  of  29  June  1980  (Figures 
5  and  6)  may  therefore  be  explained  by  large  local  scattering  but  small  regional  Rg-to-P 
scattering  because  of  the  explosion  being  fairly  deep  and  in  the  SW  Shagan  region.  The  40 
values  of  source  term  P/P-coda.  indicated  for  each  txplosion  in  Figure  1,  suggest  that  the 
low-frequency  nulls  on  the  spectra  of  most  shots  are  associated  with  relatively  low  values  of 
P/P-coda  (see  especially  spectra  of  shot  nos.  3  and  >6  in  Figure  la.  15  and  9  in  Figure  lb, 
and  7,  38,  2,  and  4  in  Figure  1c).  The  probable  e>  planation  is  that,  in  most  cases,  factors 
such  as  shallower  shot  depth  and  lower  medium  velocity  that  enhance  local  scattering  (which 
is  responsible  for  the  low-frequcncy  nulls)  also  enhance  the  generation  of  Rg  and  its  contribu¬ 
tion  by  scattering  to  P  coda. 

3.2.  ANALYSIS  OF  NORESS  ARRAY  DATA 

Analysis  similar  to  that  for  the  GKA  array  data  was  carried  out  for  the  short-period, 
vertical-component  data  from  20  Shagan  River  explosions  (Table  1)  well  recorded  at  the 
NORESS  array  which  has  25  vertical  component  sensors  located  within  a  3  km  diameter 
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aperture.  The  epicentral  distance  of  Shagan  River  explosions  to  NORESS  is  about  38°  along 
the  azimuth  of  about  313°  or  N47°W.  Again,  the  use  of  multichannel  data  and  a  least  squares 
inversion  method  isolated  the  source  and  receiver  terms  for  each  frequency.  The  derived 
source  terms  appeared  to  show  characteristics  similar  to  those  observed  for  the  EKA  data.  By 
using  P  and  contiguous  P-coda  windows  of  6.4  and  12.8  sec,  respectively  and  applying  Pazen 
taper,  spectral  ratios  P/P-coda  were  obtained  for  data  from  each  sensor  of  the  NORESS  array. 
These  ratios  were  also  used  in  a  least  squares  inversion  procedure  to  obtain  the  source  terms. 
Results  fo  two  large  (mb  =  6.0)  explosions  are  shown  in  Figure  9.  Both  figures  indicate 
spectral  ni  11s  at  the  expected  frequency  of  about  4  Hz  in  addition  to  prominent  nulls  at  fre¬ 
quency  of  about  3  Hz.  The  remarkable  similarity  of  spectral  nulls  between  spectra  of  P(6.4 
sec)  and  1  (6.4  sec)/P-coda  (12.8  sec)  for  frequency  greater  than  about  2  Hz  shows  that  the 
spectral  n  ills  in  P  are  almost  completely  filled  in  P  coda  and  P/P-coda  may  provide  an  alter¬ 
nate  effective  method  for  detecting  spectral  nulls  due  to  cancellation  by  pP. 

3.3,  ANALYSIS  OF  MllLTI-ST ATION  DATA  FOR  USSR  JVE  SHOT 

The  USSR  Joint  Verification  Experiment  (JVE)  shot  of  14  September  1988  in  SW 
Shagan  test  region  was  well  recorded  at  many  stations.  Data  from  1 1  teleseismic  (A  >  30°) 
stations  (Table  2)  was  used  to  derive  the  average  spectra  of  P(6.4  sec)  and  average  spectral 
ratio  P(6.4  sec)/P-coda  (12.8  sec).  The  average  spectra  (not  corrected  for  instrument 
response),  scaled  to  the  average  amplitudes  in  the  frequency  range  of  0.5  to  3.0  Hz,  and  the 
1 1 -station  average  ratio  o  P/P-coda,  with  no  se  lling,  are  shown  in  Figures  10  a  and  b,  respec¬ 
tively.  No.e  that  the  average  spectra  of  both  P  and  P/P-coda  indicate  nulls  at  the  expected 
frequency  >f  about  4  Hz  and  there  is  not  much  evidence  for  a  spectral  null  at  about  3  Hz 
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Figure  9.  Source  tenns  for  (a,c)  P  and  (b,d)  P/P-coda  for  two  shots  derived  by  least  squares 
inversion  of  NORESS  array  data  for  20  Shagan  River  shots.  Note  prominent  spectral  nulls  at 
about  4  He  (denoted  by  arrows)  and  3  Hz  in  both  P  and  P/P-coda  spectra  of  both  shots. 
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Figure  10.  Spectra  of  (a)  P  and  (b)  P/P-coda  for  the  Shagan  River  shot  of  14  September  1988 
derived  from  teleseismic  data  from  13  stations.  Null  at  the  expected  frequency  of  about  4  Hz 
is  not  distinct  in  the  spectra  of  P  but  is  somewhat  better  seen  in  the  spectra  of  P/P-coda. 
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observed  in  the  spectral  data  from  many  stations  such  as  the  EKA  array  and  NORESS. 


TABLE  2 

STATIONS  RECORDING  USSR  JVE  SHOT  OF  14  SEPTEMBER  1988 


No. 

STATION 

LAT 

LON 

DISTANCE 

(DEG) 

AZIMUTH 

(DEG) 

1 

ANTO 

39.900  N 

32.783  E 

33.6 

271 

2 

CHTO 

18.790  N 

98.977  E 

35.0 

145 

3 

KONO 

59.649  N 

9.598  E 

39.3 

311 

4 

GRFO 

49.692  N 

11.205  E 

42.3 

297 

5 

TOL 

39.881  N 

4.049  W 

56.7 

294 

6 

LEM 

6.833  S 

107.617  E 

61.8 

147 

7 

LON 

46.750  N 

121.810  W 

82.1 

14 

8 

SCP 

40.795  N 

77.865  W 

87.4 

342 

9 

NWAO 

32.927  S 

117.237  E 

89.2 

148 

10 

CMB 

38.035  N 

120.385  W 

90.8 

15 

11 

ANMO 

34.946  N 

106.457  W 

95.4 

4 

3.4.  ANALYSIS  OF  REGIONAL  DATA  AT  WMQ 

We  examined  the  spectra  of  regional  phases  from  17  Shagan  River  underground  nuclear 
explosions  recorded  on  the  broadband  instrument  at  the  Chinese  network  station  Uru  nchi 


(WMQ),  located  about  8.6°  or  950  km  southeast  of  the  test  site.  Tto  explosions  are  listed  in 
Table  1,  which  includes  data  from  Ringdal  and  Marshall  (1989)  for  14  Shagan  River  explo¬ 


sions  prior  to  the  year  1989.  Examples  of  seismograms  are  shown  n  Gupta  et  al.  (1990a). 
The  first  arrival  Pn  was  distinct  on  all  records,  and  the  beginning  of  Lg  was  assumed  to  be  at 


a  group  velocity  of  3.5  km/sec,  in  agreement  with  the  observed  data.  Spectra  of  Pn  were 
obtained  for  both  P(6.4  sec)  and  P(12.8  sec)  windows  and  b  /  appiying  both  Parzen  and  a  10% 
cosine  taper.  The  resulting  spectra,  without  any  smoothing,  are  shown  in  Figures  11  and  12 
for  explosions  of  13  February  1988  and  14  September  1988  (JVE),  respectively.  An 
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(a)  Pn  (6.4  sec,  Parzen  Taper)  (b)  Pn  (6.4  sec,  10%  Cosine  Taper) 


(c)  Pn  (12.8  s  c,  Parzen  Taper)  (d)  Pn  (  12.8  sec.  10%  Cosine  Taper) 


000  200  4  00  6  0C  800  1000  0  00  2  00  4  00  6  00  SOC  1000 


FREQUENCY  Hz  FREQUENCY  Hz 

Figure  11.  Spectra  and  waveforms  of  Pn  for  two  different  window  lengths  and  two  tapers  for 
the  SW  Shagan  explosion  of  13  February  1988  (mb  =  6.0)  recorded  at  WMQ.  Prominent 
nulls  in  spectra  at  about  4  Hz  are  denoti  d  by  arrows. 
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Figure  12.  Similar  to  Figure  11  for  the  explosion  of  14  September  1988  (mb  =  6.0). 
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examination  of  the  eight  spectra  of  the  two  large  explosions  confirms  the  presence  of  spectral 
null  >  at  about  4  Hz,  in  agreement  with  the  expected  value.  Note  that  Parzen  taper  provides 
better  resolution  of  the  spectral  nulls  than  the  cosine  taper.  The  same  is  true  of  shorter  versus 
longer  windows,  perhaps  because  of  less  contamination  by  later  (secondary)  arrivals  which 
will  tend  to  fill  in  the  spectral  holes. 

Using  a  window  length  of  51.2  sec  for  Lg,  the  spectra  of  Lg  and  the  spectiral  ratios 
Pn/Lg  for  explosions  of  13  February  1988  and  14  September  1988  are  shown  in  Figures  13 
and  14,  respectively.  The  spectral  ratios  are  corrected  for  noise  by  using  a  sample  of  noise 
before  the  onset  of  Pn.  In  order  to  accentuate  the  spectral  nulls  in  Pn,  the  spectral  ratios 
Pn/Lg  were  derived  by  using  only  those  data  points  for  which  the  power  ratio  of  S/N  >  1.0 
for  Pn  and  >  1.5  for  Lg.  As  expected,  the  shorter  Pn  window  and  the  use  of  Parzen  taper 
(Figures  13b  and  14b),  provided  the  most  credible  evidence  for  nulls  due  to  cancellation  by 
pP  at  the  expected  frequency  of  about  4  Hz. 

Several  Shagan  River  explosions  were  recorded  at  both  NORESS  and  WMQ.  The  spec¬ 
tra  of  teleseismic  P  and  regional  Pn  for  two  closely-spaced  explosions  are  shown  in  Figures 
15  and  16.  For  each  station,  the  spectral  shapes  including  peaks  and  nulls  of  the  two  explo¬ 
sions  appear  similar  for  frequencies  lower  than  about  3  Hz.  A  possible  reason  is  that  the 
lower  frequency  nulls  are  caused  mainly  by  near-source  scattering  which  should  be  similar  for 
closely-spaced  shots,  'fhe  NORESS  spectra  for  the  sensor  NRAO  show  a  prominent  null  at 
about  4  Hz  for  the  larg  r  explosion  of  20  June  1987  of  mb  =  6.0  and  nulls  at  about  3  Hz  for 
both  shots.  The  specu  i  of  Pn  at  WMQ  also  indicate  a  distinct  null  at  about  4  Hz  foi  the 
larger  explosion  but  thme  is  not  much  evidence  for  nulls  at  about  3  Hz.  The  Pn  spectra  of 
the  shots  of  13  February  1988  and  14  September  1988  recorded  at  WMQ  (Figures  11  and  12) 
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Parzen  Taper 


Figure  14.  Similar  to  Figure  13  for  die  SW  Shagan  explosion  of  14  September  1988. 
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6740.9  NM  5.11  SEC  3515.7  NM 

Figure  15.  Spectra  and  waveforms  of  P  at  NORESS  and  Pn  at  WMQ  for  the  SW  Shagan 
explosion  of  20  June  1987  (mb  =  6.0).  Prominent  nulls  in  spectra  at  about  4  Hz  are  denoted 
by  arrows. 


NORESS  P  (12.8  SEC,  PARZEN  TAPER)  -  ,  (b)  WMQ  Pn  (12.8  SEC.  PARZEN  TAPER) 


waveforms  of  P  at  NORESS 


also  did  not  indicate  distinct  nulls  at  about  3  Hz.  On  the  other  hand,  most  shots  recorded  at 

> »  •  '  r  *  '  i  .  *  -  *  -v  t  , .  .  •  -  *■  *  *  *  ‘  r  >  , 

EKA  (Figure  1)  showed  prominent  nulls  at  about  3  Hz.  It  seems  therefore  that  the  3  Hz  null 
is  due  to  teleseismic  path  effect,  such  as  a  multiple  arrival.  Additional  support  for  this  is  pro¬ 
vided  by  the  two  inter-shot  spectral  ratios  in  Figure  17.  The  common  3  Hz  null  in  the 
NORESS  spectra  of  both  shots  has  disappeared  whereas  the  4  Hz  null  due  to  cancellation  by 
pP  in  the  spectra  of  the  explosion  of  20  June  1987  is  prominent  because  the  pP  null  for  the 
smaller  explosions  occurs  at  a  frequency  significantly  different  than  4  Hz.  The  inter-shot 
spectral  ratio  from  WMQ  (Figure  17b)  also  provides  clear  evidence  of  a  null  at  about  4  Hz 
for  the  larger  explosion.  Note  that  the  null  observed  at  WMQ  occurs  at  a  frequency  slightly 
higher  that  of  the  null  observed  at  NORESS.  This  agrees  with  what  is  theoretically  expected 
because  of  the  somewhat  different  take-off  angles  for  Pn  and  teleseismic  P. 

3.5.  SCALING  OF  P-WAVE  SPECTRA 

It  is  important  to  study  how  the  P-wave  spectra  of  Shagan  River  explosions  vary  with 
explosive  yield  or  mb  and  to  determine  if  their  scaling  differs  significantly  from  the  scaling 
relationships  of  nuclear  explosions  proposed  by  Mueller  and  Murphy  (1971)  and  von  Seggem 
and  Blandford  (1972).  A  detailed  study  of  a  large  number  of  Shagan  River  explosions  by 
Stewart  (1988)  found  that  the  observed  P  spectra  varied  more  slowly  with  m5  than  predicted 
by  either  of  these  two  scaling  relationships  and  the  differences  are  mere  pronounced  at  lower 
frequencies  than  at  the  higher  frequencies. 

3.5.1.  Comparison  of  Theory  and  Observed  Regional  Data 

In  order  to  compare  the  observed  dependence  of  P-wave  spectra  on  yield  or  mb  with  that 
predicted  by  von  Seggem  and  Blandford  (1972)  and  Mueller  and  Murphy  (1971)  scaling 
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Figure  17.  Inter-shot  spectral  ratio  (20  June  1987)/(12  March  1987)  for  (a)  P  at  NORESS, 
and  (b)  Pn  at  WMQ  derived  from  the  spectra  in  Figures  15  and  16.  Prominent  nulls  at  about 
4  Hz  are  denoted  by  arrows. 
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relationships,  synthetic  seismograms  were  made  for  several  assumed  t*  and  pP  reflection 
coefficient  values.  As  discussed  by  Shumway  and  Blandford  (1978),  a  realistic  approach  is  to 
consider  the  pP  reflection  coefficient,  r,  to  be  frequency-dependent  and  given  by 

r  (f)  =  -  0,5  (1  +  e“f  * )  (8) 

This  relationship  implies  r(f)  =  -1  and  -0.684  at  frequencies  of  0  and  1  Hz,  respectively  and 

r(f)  approaches  -0.5  for  high  frequencies.  The  Mueller  and  Murphy  (1971)  scaling  relation¬ 
ship  depends  on  both  shot  depth  and  shot  medium  whereas  the  von  Seggern  and  Blandford 
(1972)  relationship  is  independent  of  the  shot  depth.  Shot  depths  for  a  given  value  of  mb 
were  estimated  by  using  equation  (7).  Assuming  the  P-wave  velocity  for  the  uppermost  0.5 
km  of  the  crust  in  the  Shagan  River  region  to  be  5  km/sec,  pP  and  P  are  separated  in  time  by 


tj  -■  2  sec 
a  5000 


(9) 


Synthetic  seismograms  were  constructed  for  mb  =  4.5,  5.0,  5.5,  6.0,  and  6.5  and  assuming  the 
explosion  source  medium  to  be  granite.  Three  combinations  of  t*  and  t(pP  *  P)  used  were: 


(1)  t*  =0.1.  t(pP-P)  =  td, 

(2)  t*  =  0.1,  t(pP  -  P)  =  2  td,  and 

(3)  t*  =  0.2,  t(pP  -  P)  =  2  td. 


The  use  of  t*  =  0.1  for  WMQ  may  be  a  good  approximation  on  the  basis  of  Sereno’s  (1990) 
study  suggesting  a  frequency-independent  Q  value  of  1175.  The  instrument  response  was  that 
of  the  broadband  short-period  instrument  at  WMQ. 

Synthetic  seismograms  based  on  the  use  of  both  von  Seggern  and  Blandford  (1972)  and 
Mueller  and  Murphy  (1971)  scaling  relationship  and  the  observed  data  were  bandpassed  for 
the  frequency  bands  of  0.5- 1.0  Hz,  4.0-6.0  Hz,  and  7.0-9.0  Hz  by  using  a  three-pole  phaseless 


Butterworth  filter.  For  the  observed  data,  only  14  Shagan  River  explosions  prior  to  the  year 
1989  for  which  precise  mb  values  were  available  from  Ringdal  and  Marshall  (1989)  were 
used.  Log  RMS  values  were  computid  for  the  P(12.8  sec)  windows  with  10%  cosine  taper 
and  correction  for  instrument  respon  e.  For  the  observed  data,  a  correction  for  noise  was 
made  by  using  a  sample  of  noise  prior  to  the  onset  of  P  and  subtracting  the  power  in  noise 
from  that  in  signal.  Plots  of  log  RMS  values  versus  mb  for  the  synthetic  and  the  observed 
data  for  frequency  b;  nds  of  0.5- 1.0  Hz,  4.0-6.0,  and  7.0-9.0  are  shown  in  Figures  18,  19,  and 
20,  respectively.  Th<  amplitude  scale  is  arbitrary.  The  three  sets  of  synthetics  provide  results 
that  are  not  much  different  whether  von  Seggem  and  Blandford  or  Mueller  and  Murphy  model 
is  used.  There  is,  however,  a  large  discrepancy  between  the  observed  and  synthetic  data  with 
the  observations  indicating  much  slower  variation  with  mb  for  the  lowest  frequency  band  of 
0.5- 1.0  Hz  (Figure  18).  The  discrepancy  is  considerably  reduced  for  the  frequency  band  of 
4.0-6.0  Hz  (Figure  19)  and  appears  further  reduced  for  the  frequency  band  of  7.0-9.0  Hz  (Fig¬ 
ure  20). 

The  results  for  spectral  ratios  (0.5-1.0  Hz)/(4.0-6.0  Hz)  and  (0.5-1.0  Hz)/(7.0-9.0  Hz)  for 
the  two  theoretical  modeh  and  the  observed  data  are  shown  in  Figures  21,  22,  23,  and  24  in 
which  the  amplitude  scale  is  not  arbitrary.  A  comparison  of  the  observed  and  the  theoretical 
results  again  indicates  the  low-frequency  spectra  to  vary  more  slowly  with  mb  than  expected 
from  either  of  the  two  theoretical  models. 

3.5.2.  Comparison  of  Theory  and  Observed  Teleseismic  Data 

A  comparison  of  the  scaling  of  spectra  with  r.  .b  between  the  theoretical  results  based  on 
von  Seggem  and  Blandford  (1972)  and  Mueller  and  Murphy  (1971)  and  observation  was 
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Figure  18.  Log  RMS  versus  mb  for  frequency  band  of  0.5- 1.0  Hz  for  observed  data  at  WMQ 
and  three  synthetics  with  t*  and  t(pP-P)  values  as  indicated.  The  amplitude  scale  is  arbitrary. 
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(b)  MUELLER  AND  MURPHY  (1971) 
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(a)  VON  SEGGERN  AND  BLANDFORD  (1972) 
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Figure  19  Log  RMS  versus  m*,  for  frequency  band  of  4.0-6.0  Hz  for  observed  data  at  WMQ 
and  three  ynthetics  with  t*  and  t(pP-P)  values  as  indicated.  The  amplitude  scale  is  arbitrary. 
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(a)  VON  SEGGERN  AND  BLANDFORD  (1972)  (b)  MUELLER  AND  MURPHY  (1971) 


A  t*  =  0.1,t(pP-P)  =  td 
A  t‘ =  0.1,t(pP-P)  =  2td 
0  t*  =  0.2,  t(pP  -  P)  =  2  td 

Figure  20.  Log  RMS  versus  ntt,  for  frequency  band  of  7.0-9.0  Hz  for  observed  data  at  WMQ 
and  three  synthetics  with  t*  and  t(pP-P)  values  as  indicated.  Tht  amplitude  scale  is  arbitrary. 
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Figure  21.  Log  RMS  ratio  (0.5-1.0)/(4.0*6.0)  versus  mb  for  observed  data  at  WMQ  and  three 
synthetics  based  on  von  Seggem  and  Blandford  (1972)  scaling  and  t*  and  t(pP-P)  values  as 
indicated. 
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Figure  22.  Log  RMS  ratio  (0.5-1.0)/(4.0-6.0)  versus  mb  for  observed  dat;  at  WMQ  and  three 
synthetics  based  on  Mueller  and  Murphy  (1971)  scaling  and  t*  and  t(pP*P)  values  as  indi¬ 
cated. 
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Figure  23.  Log  RMS  ratio  (0.5-1.0)/(7.  )-9.0)  versus  mj,  for  observed  data  at  WViQ  and  three 
synthetics  based  on  von  Seggern  and  I  iandford  (1972)  scaling  and  t*  and  t(pP-P)  values  as 
indicated. 
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Figure  24.  Log  RMS  ratio  (0.5*  1.0)/(7. 0*9.0)  versus  mb  for  observed  data  at  WMQ  and  three 
synthetics  based  on  Mueller  and  Murphy  (1971)  scaling  and  t*  and  t(pP-P)  values  as  indi¬ 
cated. 
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made  by  using  teleseismic  data  from  25  Shagan  River  explosions  recorded  at  the  center  ele¬ 
ment  of  NORSAR  subarray  1A.  Most  of  these  data  had  good  S/N  up  to  about  J  Hz.  Syn¬ 
thetic  seismograms  were  again  made  for  the  same  two  assumed  values  of  t*  and  the  same  two 
values  of  t(pP  -  P).  In  fact,  all  parameters  were  selected  to  be  the  same  except  that  a  Parzcn 
taper  was  used  on  6.4  sec  long  P-windows  and,  of  course,  NAO  instrument  response  was 
used.  Note  that  Der  etal.  (1985)  suggested  t*  =  0.14  for  Shagan  River  explosions  recorded 
at  NORSAR.  The  results  for  spectral  ratios  (0.5-2.0  Hz)/(3.0-5.0  Hz)  and  (0.5-2.0  Hz)/(7.0- 
9.0  Hz)  for  the  two  theoretical  models  and  the  observed  data  are  shown  in  Figures  25,  26,  27, 
and  28.  A  companion  of  the  observed  and  the  theoretical  results  shows  large  discrepancy, 
especially  for  the  smaller  mb  shots.  The  discrepancy  is  of  the  same  type  as  for  the  regional 
data  (Figures  21  through  24)  but  somewhat  larger  in  magnitude.  These  results  also  suggest 
that  the  observed  low-frequency  spectra  vary  more  slowly  with  mb  than  expected  from  either 
of  the  two  theoretical  models. 

3.6.  DISCUSSION 

A  three-dimensional  geological  model  of  the  Shagan  River  test  site  by  Leith  and  Unger 
(1989)  shows  a  distinct  difference  betw  een  the  NE  and  SW  portions  of  the  test  site,  with  the 
granites  closer  to  the  surface  and  the  alluvium  thinner  in  the  southwest.  Finite  difference 
simulations  of  scattering  due  to  ne;  source  geology  for  the  Yucca  Flats  explosions  by 
McLaughlin  etal.  (1987)  ami  Stead  ;  nd  Helmberger  (1988)  indicate  large  amplitude  arrivals 
within  a  few  sec  of  the  first  I*.  P-wavc,  spectra  of  the  Pahute  Mesa  explosion,  MAST,  shows 
a  spectral  null  at  the  frequency  expected  for  interference  by  pP  for  short  time  windows  and  at 
considerat  y  lowet  frequency  for  longer  windows  (Gupt;  etal.,  1991).  F-k  analysis  of 
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Figure  25.  Log  RMS  ratio  (0.5-2.0V(3.0-5.0)  versus  mb  for  observed  data  at  NORSAR  and 
three  synthetics  based  on  von  Seggem  and  Blandfotd  (1972)  scaling  and  t*  and  t(pP*P)  values 
as  indicated. 
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Figure  26.  Log  RMS  ratio  (0.5-2.0)/(3.0-5.0)  versus  mb  for  observed  data  at  NORSAR  and 
three  synthetics  based  on  Muelier  a  id  Murphy  (1971)  scaling  and  t*  and  t(pP-P)  values  as 
indicated. 
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Figure  27.  Log  RMS  ratio  (0.5-2.0y(7.0-9.0)  versus  mb  for  observed  data  at  NORSAR  and 
three  synthetics  based  on  von  Seggem  and  Blandford  « 1972)  scaling  and  t*  and  t(pP-P)  values 
as  indicated. 
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Figure  28.  Log  RMS  ratio  (0.S-2.0)/(7. 0-9.0)  versus  for  observed  data  at  WMQ  and  three 
synthetics  based  on  Mueller  aid  Murphy  (1971)  scaling  and  t*  and  t(pP-P)  values  as  indi¬ 
cated. 
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teleseismic  records  of  Yucca  Flat  explosions  by  Gupta  etal.  (1990b)  suggested  Rg-to-P 
arrivals  within  a  few  sec  of  the  direct  P.  A  finite  difference  study  of  scattering  due  to  irregu¬ 
lar  low-velocity  surface  layers  by  Levander  and  Hill  (1985)  showed  that  when  a  planar  pulse 
is  incident  upon  the  low-velocity  layer,  the  roughness  induces  a  strong  resonant  coupling  to 
Rayleigh  modes,  a  condition  which  cannot  occur  for  a  laterally  invariant  medium.  The 
tesonant  coupling  causes  a  long  coda  following  the  direct  seismic  arrival.  Such  a  mechanism 
may  be  responsible  for  the  scattering  of  explosion-generated  Rg  into  teleseismic  P  coda. 
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4.  STUDY  OF  NOVAYA  ZEMLYA  EXPLOSIONS 


As  for  most  other  test  sites,  there  is  considerable  controversy  over  the  interpretation  of 
pP  parameter  estimates  of  Novaya  Zemlya  explosions  (e.  g.  Burdick,  1990;  Lay,  1991). 
There  is  general  agreement  over  the  existence  of  strong  pP  reflections  but  considerable 
disagreement  regarding  the  pP-P  delay  times  as  well  as  pP/P  amplitudes.  The  most  reliable 
delay  times,  consistent  with  the  geology  of  the  test  site,  have  so  far  been  those  derived  by  the 
maximum  likelihood  mutichannel  deconvolution  methods  (Chan  etal.,  1988;  Der  etal., 
1987b).  We  used  inter-shot  spectral  ratios,  averaged  over  all  available  sensors,  to  obtain 
significantly  more  precise  delay  times  for  several  large  shots.  Burdick  (1990)  showed  that  the 
spectral  estimation  procedure  used  in  the  spectral  averaging  method  of  Murphy  (1989)  tends 
to  underestimate  the  pP  amplitudes  due  to  spectral  smoothing.  We  made  a  comparison  of 
theoretical  and  observed  source  spectra  in  order  to  estimate  pP/P  amplitudes  of  several  large 
explosions. 

4.1.  P-WAVE  SOURCE  SPECTRA  FROM  U.K.  ARRAY  DATA 

The  U.  K.  array  data  from  the  EKA  (A  =  29°,  source  azimuth,  d>  =  264°),  GBA 
(A  =  61°,  <I>  =  155°),  and  YKA  (A  =  41°,  <I>  =  353°)  were  collected  for  as  many  as  19  Novaya 
Zemlya  explosions  (Table  3)  out  of  which  only  one  explosion  (27  September  1973  or  73270) 
is  from  the  southern  Novaya  Zemlya  test  site  region.  The  number  of  usable  (free  from  clip¬ 
ping  or  spike)  sensors  available  at  each  array  for  each  explosion  is  indicated  in  the  table.  For 
the  largest  explosion  (66300),  all  channels  of  data,  except  one,  were  clipped.  Most  of  the  data 
in  Table  3  is  from  Lilwall  and  Marshall  (1986)  who  carried  out  a  least  squares  joint  epicenter 
estimate  of  origin  time  and  epicenter  together  with  a  maximum-likelihood  estimate  of 
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magnitude.  The  surface  topography  in  the  northern  Novaya  Zemlya  region  is  characterized  by 
steep  mountain  slopes  whereas  topography  at  the  southern  site  is  less  pronounced. 


TABLE  3 

NOVAYA  ZEMLYA  SHOTS  RECORDED  AT  UK  ARRAYS 


No. 

DATE 

LAT 

LON 

ORIGIN  TIME 

mb 

NUMBER 

EKA 

OF  SENSORS 

GBA  YKA 

1 

25  Oct  1964 

64299 

73.386 

54.997 

07:59:58.1 

4.82 

18 

0 

0 

2 

27  Oct  1966 

66300 

73.388 

54.845 

05:57:58.1 

6.47 

1 

0 

0 

3 

07  Nov  1968 

68312 

73.388 

54.873 

10:02:05.5 

6.11 

1 

1 

2 

4 

14  Oct  1969 

69287 

73.389 

54.796 

07:00:06.6 

6.18 

2 

1 

0 

5 

27  Sep  1973 

73270 

70.756 

53.746 

06:59:58.5 

5.83 

16 

0 

17 

6 

29  Sep  1976 

76273 

73.360 

54.880 

02:59:57.7 

5.77 

1 

2 

17 

7 

20  Oct  1976 

76294 

73.399 

54.835 

07:59:58.1 

4.89 

20 

17 

0 

8 

01  Sep  1977 

77244 

73.339 

54.626 

02:59:58.0 

5.71 

8 

0 

17 

9 

09  Oct  1977 

77282 

73.414 

54.935 

10:59:58.1 

4.51 

15 

9 

0 

10 

10  Aug  1978 

78222 

73.293 

54.885 

07:59:58.0 

6.04 

0 

0 

10 

11 

27  Sep  1978 

78270 

73.350 

54.677 

02:04:58.6 

5.68 

15 

1 

18 

12 

24  Sep  1979 

79267 

73.346 

54.679 

03:29:58.8 

5.80 

19 

0 

17 

13 

18  Oct  1979 

79291 

73.318 

54.821 

07:09:58.8 

5.85 

20 

1 

17 

14 

11  Oct  1980 

80285 

73.335 

54.938 

07:09:57.5 

5.80 

0 

10 

14 

15 

01  Oct  1981 

81274 

73.308 

54.817 

12:14:57.3 

5.91 

0 

11 

11 

16 

11  Oct  1982 

82284 

73.348 

54.601 

07:14:58.7 

5.52 

15 

16 

17 

17 

18  Aug  1983 

83230 

73.358 

54.974 

16:09:58.9 

5.84 

0 

19 

16 

18 

25  Sep  1983 

83268 

73.326 

54.564 

13:09:58.2 

5.71 

12 

0 

17 

19* 

25  Oct  1984 

84299 

73.370 

54.960 

06:29:57.7 

5.90 

0 

16 

13 

*  from  Sykes  and  Ruggi  (1986) 
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Frequency-domain  source  terms  were  derived  by  a  least  squares  analysis  of  all  available 
data  from  each  array  for  the  Novaya  Zemlya  explosions  listed  in  Table  3.  Results  from  14 
explosions  recorded  at  the  EKA  array,  based  on  use  of  the  P(6.4  sec)  window,  are  shown  in 
Figure  29  in  which  the  source  spectra  are  arranged  in  order  of  decreasing  mb.  Note  that 
source  spectra  for  the  three  largest  explosions  are  based  on  extremely  limited  data  (see  Table 
3)  so  that  the  corresponding  source  spectra  should  be  less  reliable.  Results  of  similar  analysis 
carried  out  by  using  the  P(12.8  sec)  windows  are  shown  in  Figure  30.  An  examination  of  the 
source  spectra  in  Figures  29  and  3(  indicates  rather  distinct  spectral  nulls  (indicated  by 
arrows)  for  most  of  the  larger  shots  at  a  frequency  of  about  2.5  to  3.0  Hz.  The  spectral  nulls 
in  I  igure  30  are  generally  less  distinct  than  in  Figure  29,  probably  due  to  greater  contamina¬ 
tion  by  later  arrivals  in  the  1(  nger  P  window.  Seismic  velocities  in  the  shallow  crust  at  the 
Novaya  Zemlya  test  site  are  somewhat  smaller  than  those  at  the  Shagan  River  test  site  (V. 
Adushkin,  oral  communication).  Since  the  spectral  nulls  due  to  pP  for  large  Shagan  River 
explosions  occur  at  a  frequency  of  about  4  Hz,  the  spectral  nulls  observed  at  frequency  of  2.5 
to  3.0  Hz  for  the  larger  Novaya  Zemlya  shots  are  also  expected  to  be  due  to  cancellation  by 
pP. 

The  GBA  data  for  12  Novaya  Zemlya  explosions  (Table  3)  was  analyzed  in  the  same 
manner  as  the  EKA  data  and  the  soui  :e  spectra,  based  on  use  of  the  P(6.4  sec)  window,  are 
shown  in  Figure  3 1 .  The  spectra  are  arranged  in  order  of  decreasing  mb.  A  few  of  the  larger 
shots  show  prominent  spectral  nulls  (indicated  by  arrows)  at  frequency  that  lies  between  2.5 
and  3.0  Hz.  Results  of  similar  analysis  of  the  YKA  data  from  14  explosions  are  shown  in 
Figure  32.  Again,  the  source  spectra  of  many  larger  shots  suggest  prominent  spectral  nulls  at 
frequency  of  about  2.5  to  3.0  Hz.  It  should  be  noted  that  the  spectral  nulls  for  a  given 
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Figure  29.  Source  spectra  (not  corrected  for  instrument  response)  of  14  Novaya  Zemlya 
explosions  derived  from  P(6.4  sec)  window  of  EKA  array  data.  The  shots  are  arranged  in 
order  of  decreasing  mb.  Prominent  spectral  nulls  are  indicated  by  arrows. 
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Figure  30.  Similar  to  Figure  29  for  the  P(12.8  sec)  window.  The  spectral  nulls  are  generally 
less  distinct  than  in  Figure  29  due  to  contamination  by  later  arrivals. 
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Figure  31.  Similar  to  Figure  29  for  12  explosions  recorded  at  the  GBA  array. 
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Figure  32.  Similar  to  Figure  29  for  14  explosions  recorded  at  the  YKA  array. 
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explosion  observed  at  the  three  arrays  (Figures  29,  31,  and  32)  are  generally  at  somewhat 
different  frequencies.  Possible  reasons  are  different  -.s  in  source-receiver  paths  and  take-off 
polar  and  azimuthal  angles. 

Spectral  ratios  P/P-coda,  with  a  P  window  of  6.4  sec  and  the  contiguous  P-coda  window 
of  12.8  sec,  were  computed  for  all  usable  data  from  th  e  three  arrays  (Table  1).  As  mentioned 
earlier,  the  pP  nulls  in  spectra  of  P-coda  are  expected  to  be  significantly  less  distinct  than  in 
spectra  of  the  initial  P  so  that  the  spectral  ratio  P/P-coda  should  be  useful  for  determining  the 
null  frequencies.  In  order  to  emphasize  the  spectral  mills,  power  thresholds  of  S/N  >  1.0  and 
l.S  were  applied  to  the  P  and  P-coda  spectra,  respectively.  Using  a  least  squares  regression 
technique  applied  to  each  spectral  amplitude,  frequency-domain  source  terms  were  obtained 
for  all  explosions  except  shot  66300  for  v  hich  data  from  only  one  sensor  is  available  (Table 
3).  The  source  terms  P/P-coda  for  18  explosions,  shown  in  Fi  gure  33,  indicate  distinct  spec¬ 
tral  nulls  (indicated  by  arrows)  for  most  shots  of  mb  larger  than  5.5.  Null  frequencies  for  the 
large  shots  are  at  about  2.5  to  3.0  Hz  and  do  not  show  a  systematic  variation  with  mb. 

4.2.  INTER-SHOT  RATIOS  AND  ESTIMATES  OF  NULL  FREQUENCY 

Spectral  ratios  of  closely  spaced  explosions  at  a  common  recording  station  provide 
source  information  that  is  virtually  free  from  propagation  path  effects,  including  attenuation 
and  receiver  function  (King  ei  al. ,  1972).  The  procedure  can  be  used  to  estimate  spectral 
nulls  or  delay  times  if  shot  depths  of  the  two  explosions  are  significantly  different.  Numerous 
applications  of  this  technique  have  indicated  the  pP  delay  time  to  be  longer  than  expected 
from  known  overburden  velocities  and  its  amplitude  significantly  smaller  than  th  it  expected 
for  elastic  free  surface  reflection. 
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Figure  33.  (Continued). 


Inter-shot  spectral  ratios  were  obtained  by  combining  a  large  shut  with  a  small  shot 
recorded  at  common  stations.  P(6.4  sec)  long  windows  were  used  and,  in  order  to  emphasize 
nulls  in  spectra  of  the  larger  shot,  a  power  ratio  of  S/N  >  1.0  for  the  larger  explosion  and  S/N 
>  1.5  for  the  smaller  explosion  were  used.  Results  from  the  EKA  data  for  the  four  largest 
shots  (each  with  3  or  more  usable  sensors),  combining  each  with  three  smallest  available 
shots,  are  shown  in  Figures  34  and  35.  These  plots  are  based  on  the  average  of  all  common 
sensors  and  only  those  frequency  values  for  which  data  from  three  or  more  sensors  were 
available  are  shown.  The  vertical  lines  represent  one  standard  deviation  error  bars  derived 
from  the  values  from  3  or  more  sensors.  Distinct  spectral  nulls  are  observed  on  all  plots  in 
Figures  34  and  35  and  the  null  frequency  for  each  of  the  four  large  explosions  stays  remark¬ 
ably  constant  when  the  large  shot  is  combined  with  three  different  smaller  shots.  Results 
derived  from  similar  analysis  of  data  from  two  large  explosions  recorded  at  GBA  and  YKA 
are  shown  in  Figures  36  and  37,  respectively.  In  Figures  37a,  b,  c,  error  bars  are  not  given 
since  data  from  only  two  sensors  were  available  for  shot  68312.  Distinct  and  consistent  spec¬ 
tral  nulls  arc  again  observed  at  frequencies  of  about  2.5  to  3.0  Hz  in  both  Figures  36  and  37. 

The  spectral  null  frequencies  for  8  relatively  large  explosions  in  Figures  34  through  37 
agree  well  with  those  in  Figure  33.  This  is  especially  true  of  explosions  79267  (Figures  34d, 
e,  f  and  33i),  81274  (Figures  36a.  b,  c  and  33d),  and  68312  (Figures  37a,  b,  c  and  33b).  It 
seems  therefore  that  inter-shot  spectral  ratios  provide  a  simple  and  effective  method  of 
estimating  null  frequencies  and  corresponding  delay  times. 
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(d)  79267(5.80)  /  64299(4.82) 


(c)  79291(5.85)/ 77282(4.51) 


O 


'0.0*  2.00  4.00  6.00  '0.00 


FREQUENCY  (HZ) 


(f)  79267(5.80)/ 77282(4.51) 


2.00  4^00  6.00 

FREQUENCY  (HZ) 


Figure  34.  Inter-shot  spectral  ratios,  averaged  over  all  common  sensors,  for  two  larje  explo¬ 
sions  recorded  at  EKA,  each  combined  with  three  smaller  shots.  Prominent  and  consistent 
spectral  nulls  (indicated  by  arrows)  are  observed  at  about  3  Hz.  fhe  vertical  lines  represent 
error  bars  with  one  standard  deviation. 
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Figure  35.  Similar  to  Figure  3 1  for  two  additional  larger  shots.  Distinct  and  consistent  spec¬ 
tral  nulls  are  again  observed. 
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Figure  36.  Similar  to  Figure  34  for  two  large  explosions  recorded  at  GBA  showing  distinct 
and  consistent  spectral  nulls. 
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Figure  37.  Similar  to  Figure  34  for  two  large  explosions  recorded  at  YKA  showing  consistent 
spectral  nulls  that  are  considerably  more  prominent  for  shot  68312)  than  for  shot  78222. 


4.3.  ESTIMATES  OF  pP/P  BY  COMPARISON  WITH  SYNTHETICS 


In  order  to  resolve  the  large  ambiguity  associated  with  pP/P  amplitudes  as  determined  by 
various  investigators,  we  selected  a  few  events  with  distinct  pP  spectral  nulls  and  compared 
their  spectra  with  those  of  theoretical  waveforms  (synthetics).  Three  relatively  large  explo¬ 
sions  selected  for  this  purpose  are  73270  from  southern  Novaya  Zemlya  and  79267  and  79291 
from  northern  Novaya  Zemlya.  The  reasons  for  their  selection  are:  (1)  data  from  a  large 
number  of  sensors  are  available  (Table  3),  (2)  the  EKA  source  spectr  of  these  explosions  in 
Figure  29  appear  to  provide  distinct  spectral  nulls  at  frequency  of  2.5,  2.8,  and  2.8  Hz,  respec¬ 
tively,  and  (3)  the  three  null  frequencies  agree  remarkably  well  with  those  derived  from  P/P- 
coda  (Figure  33)  and  from  inter-shot  spectral  ratios  (Figure  34).  According  to  Sykes  and 
Ruggi  (1986),  these  3  explosions  had  estimated  yields  of  100,  55,  and  70  kt,  respectively. 
Following  Chan  et  al.  (1988)  for  Novaya  Zemlya  explosions  recorded  at  EKA,  t*  was 
assumed  to  be  0.23  whereas  the  times  were  assumed  to  be  those  given  bv  the  observed  null 
frequencies.  Theoretical  seismograms  were  obtained  by  using  von  Seggern  and  Blandford 
(1972)  scaling  relationship  for  several  different  assumed  values  of  (frequency  independent) 
pP/P  amplitude.  Spectra  of  theoretical  seismt  grams  were  obtained  by  using  signal  windows 
of  P(6.4  sec)  and  P(3.2  sec)  with  Parzen  taper.  These  spectra  were  compared  with  the  source 
spectra,  derived  by  using  the  same  window  lengths,  for  the  three  explosions  mentioned  above. 
Some  of  the  results  are  shown  in  Figures  38,  39,  and  40  which  suggest  a  pP/P  value  of  about 
0.4  to  0.6  in  nearly  all  cases. 

An  examination  of  the  spectra  of  synthetics  in  Figures  38,  39,  and  40  suggests  that  the 
pP/P  amplitude  exerts  strong  influence  on  shape  of  the  spectra  at  fre  |uencies  close  to  the  pP 
null  frequency  (or  "strength"  of  the  pP  null).  Variations  in  yield  have  oniy  a  minor  influence. 
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Figure  38.  Theoretical  and  observed  spectra  for  shot  73270  recorded  at  EKA  array.  Top  two 
spectra  are  based  on  reflection  coefficien:  values  of  0.6  and  0.4,  respectively  and  the  use  of 
(a)  P(6.4  sec)  and  (b)  P(3.2  sec)  windows  Comparisons  suggest  pP/P  of  about  0.6  and  0.4  for 
the  observed  spectra  in  (a)  and  ( >),  respectively. 


67 


LOG  AMPLITUDE 

0.00  1 .00  2.00  3.00  4.00  5.00  6.00 


Figure  40.  Similar  to  I  igure  38  for  shot  79291.  Comparisons  suggest  pP/P  of  about  0.6  each 
for  the  observed  spectni  in  both  fa)  and  (b). 
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Construction  of  additional  synthetics  showed  that  reast  nable  variations  in  t*  also  did  not  have 
much  effect  on  the  strength  of  the  pP  null.  A  compar  son  of  the  P(6.4  sec)  spectra  in  Figures 
29,  31,  and  32  with  the  spectra  of  synthetics  would  herefore  suggest  pP/P  ampli  ude  to  be 
about  0.4-0.6  for  most  large  shots  and  considerably  smaller  than  0.4  for  the  smaller  mb  shots. 


70 


5.  CONCLUSIONS 


A  large  amount  of  teleseismic  and  regional  data  from  Shagan  River  explosions  have  been 
analyzed  by  using  several  methods  of  analysis  such  as  frequency-domain  isolation  of  source 
and  receiver  terms,  multichannel  deconvolution,  use  of  spectral  ratios  P/P-coda  and  Pn/Lg  for 
teleseismic  and  regional  data,  respectively,  and  imer-shot  spectral  ratios  at  common  stations. 
The  P-wave  spectra  appear  to  show  three  different  types  of  spectral  nulls:  distinct  nulls  due  to 
cancellation  by  pP  at  the  expected  frequency  of  about  4  Hz  in  both  teleseismic  and  regional 
data  for  large  shots  (mb  =  6),  low-frequency  (less  than  about  2  Hz)  nulls  in  teleseismic  data, 
generally  stronger  for  smaller  explosions  and  for  shots  in  the  NE  Shagan  region,  and  nulls  at 
about  3  Hz  in  teleseismic  data  for  most  shots.  The  low-frequency  nulls  are  probably  due  to 
near-source  scattering  whereas  the  3  Hz  null  appears  to  be  due  to  teleseismic  path  effect  such 
as  a  multiple  arrival.  An  investigation  of  the  scaling  of  P-wave  spectra  is  carried  out  by  com¬ 
parison  of  observed  data  from  WMQ  and  NORSAR  with  synthetics  assuming  several  t*  and 
pP  reflection  coefficient  values.  The  observed  P  spectra  vary  more  slowly  with  mh  than 
predicted  by  either  Mueller  and  Murphy  (1971)  or  von  Seggem  and  Blandford  (1972)  scaling 
relationships,  in  qualitative  agreement  with  Stewart’s  (1988)  analysis  of  Shagan  River  explo¬ 
sions  recorded  at  the  U.K.  arrays.  The  differences  between  theoretical  and  observed  spectra 
are  greater  at  lower  than  at  the  higher  frequencies. 

Source  terms  derived  from  P-wave  spectra  of  several  larger  Novaya  Zemlya  explosions 
show  distinct  nulls  at  frequencies  close  to  those  expected  for  nulls  due  to  pP.  The  spectral 
null  frequencies  differ  by  small  amounts  from  one  array  to  another,  presumably  because  of 
variation  in  path  effect  .  Spectral  ratio  of  P/P-coda,  averaged  over  all  available  sensors  at  all 
arrays,  provides  anothe  method  for  estii  tating  spectral  nulls.  The  simplest  and  most  reliable 
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method  for  estimating  null  frequencies  of  larger  explosions  is  found  to  be  inter-shot  spectral 
ratios  at  common  receivers.  The  frequency  of  spectral  nulls  determined  by  this  method  stays 
remarkably  constant  when  the  larger  shot  is  combined  with  various  smaller  shots  and  the 
results  agree  well  with  those  determined  by  the  spectral  ratio  P/P-coda  averaged  over  ail  sen¬ 
sors.  A  knowledge  of  pP  spectral  nulls  may  >e  used  to  determine  shot  depths  if  the  overbur¬ 
den  velocities  are  known.  Delay  times  of  large  Novaya  Zemlya  explo;  ions  determined  in  this 
study  are  in  good  agreement  with  those  from  source  deconvolutions  and  appear  to  be  con¬ 
sistent  with  the  overburden  velocities  in  the  test  site  region.  A  comparison  of  the  observed 
data  with  synthetics  suggests  the  pP/P  amplitudes  for  the  larger  events  to  be  about  0.4  to  0.6. 
The  null  frequencies  do  not  seem  to  vary  significantly  with  mb;  this  will  be  consistent  with 
Leith  et  al.  ’s  (1990)  suggestion  that  explosions  at  the  northern  test  site  are  emplaced  in  near¬ 
horizontal  tunnels  so  that  the  larger  shots  may  be  relatively  underburied  and  therefore  associ¬ 
ated  with  somewhat  larger  than  expected  pP  null  frequencies. 
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